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We report detailed and extensive first-principles molecular-dynamics (MD) simulations of the struc- 
ture and electronic properties of amorphous InP produced by rapid quenching from the liquid. The 
structure of the material is found to be strongly ordered chemically, even though there is a signifi- 
cant number of coordination defects and despite the presence of odd-membered rings. We find, as 
a consequence, that there exists "wrong bonds" in the system, in an amount of about 8%; these 
result from the presence of coordination defects, not of local composition fluctuations, as has been 
conjectured. The system, in fact, is found to be over-coordinated, which might be the reason for 
the observed higher density of a-InP compared to c-InP. We have also investigated the possibility of 
pressure-amorphizing InP. Our calculations indicate that the cost of a transformation of the com- 
pressed zinc-blende crystal into an amorphous phase is so large that it is very unlikely that it would 
take place. 

PACS numbers: 61.43.Bn, 61.43.Dq, 64.70.Pf, 71.23.Cq 
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I. INTRODUCTION 

Despite considerable work over the last few decades, 
precise understanding of the static and dynamic struc- 
ture of glasses and amorphous materials temains a chal- 
lenge to theorists and experimentalists E^B Only average 
properties of these materials are accessible to experiment; 
even in such a simple material as o-Si, a covalent semicon- 
ductor, detailed experimental knowledge of the atomic 
arrangements on the local lengthscale is missing. The 
average coordination number of a-Si, for instance, is not 
known exactly, though it appears that it is almost the 
same as for c-Si, i.e., 4.EI The only way of obtaining de- 
tailed microscopic information on the local atomic struc- 
ture is thus via theoretical modelling. In particular, ab- 
initio molecular-dynamics (MD) simulations, which de- 
scribe accurately the interatomic potentials, have been 
able to generate structural models of a-Si and a-GaAs 
which yield meas.urable quantities in good agreement 
with experiment .LfB Ab-initio calculations are, however, 
computationally very demanding. Empirical potentials 
such as Stillinger-WeberO or TersofSiJ work reasonably 
well for Si, Ge and their alloys, but there exist no such 
models for III-V compounds. Indeed, these materials 
are inherently much more difficult to model than the 
corresponding elemental systems because of the added 
complexity of (partly) ionic bonding, which results in a 
strong degree of chemical order in the crystal. Despite 
these difficulties, a set of transferable tight-binding (TB) 
potentials has been developed, for some of the III-V's, 
in partipilas GaAs and GaP.lliJ In recent MD studies of 
a-GaAst3"L3 and Z-GaAs|la these models were found to 
produce results in good agreement with those from ab- 



initio studies and from experiment. 

InP is an important material for the industry of micro- 
electronic and optoelectronic devices, in particular ipJthe 
field of high-speed computing and communications £3 In 
spite of its potential, it has been much less studied than 
GaAs (the two materials have similar bandgaps), mostly 
because of the difficulty in fabricating high-quality InP 
crystals in large enough quantities. Resorting to the 
amorphous phase of the material might be a way out 
of this problem; indeed, a-InP is expected to find its way 
in the fabrication of integrated circuits. 

Attempts at fitting a reasonable tight-binding model 
for InP have so far been unsuccessful and, as mentioned 
above, there exists no empirical or other potential for this 
material. A first-principles approach therefore seems to 
be the only possible avenue for constructing models of the 
amorphous material. Here we propose a model for sto- 
ichiometric a-InP obtained by a melt-and-quench cycle. 
To our knowledge, this constitutes the very first attempt 
at constructing a realistic model of a-InP. Experimen- 
tally, the structure of the material remains to a large 
extent unresolved, despite the fact that some structural 
measurements have been reported in the literature (see 
below). Several questions pertaining to the local atomic 
order remain open. In particular, though it is clear that 
a-InP is disordered from both chemical and structural 
viewpoints, experiment has not yet given a precise value 
for the proportion of wrong bonds in the material — and 
its relation to coordination fluctuations. Thus, for in- 
stance, it is not clear if wrong bonds result from the pres- 
ence of topological defects, such as odd-membered rings, 
or from local compositional fluctuations (i.e., antisitea or 
clustering) arising from conditions of preparation.EIl In 
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the case of InP, which is rather strongly ionic, heteropo- 
lar bonding should be favored over homopolar bonding, 
and the proportion of wrong bonds consequently reduced. 

a-InP is normally produced by flash evaporation jof jfe 
InP and deposition onto an appropriate substrate,llZrEil 
but it can also be obtained by ion implantation — which 
in principle yields better-quality material with repro- 
ducible properties — though usually not in quantity suf- 
ficient for such atooiic structural measurements as x-rays 
to be carried outHj There are only vepjj few reports of 
ion-implantation-amorphization of InP.E3 Other covalent 
semiconductors, such as Si and GaAs, cauj-alpa be pro- 
duced through proper pressure treatment .cirE3 To our 
knowledge, a-InP has never been produced in this way. 
It has however been conjectured that pressure-induced 
amorphization should pot occur in strongly-ionic com- 
pound semiconductors E3 This conjecture has not been 
verified; it is clearly of interest to examine the question 
in the case of InP, which is significantly more ionic than 
GaAs (0.421 vs 0.310 — cf. Ref. |8|). 

We have therefore also examined, in the course of this 
study, the possibility of fabricating a-InP through the ap- 
plication of pressure. Our computer simulations indicate 
that InP does not amorphize under pressure, even for val- 
ues largely in excess of those required for the system to 
transform into the high-pressure NaCl phase. The energy 
of the compressed zinc-blende crystal, we find, remains 
lower than that of the amorphous phase produced from 
the melt, even at higher density. Thus, there is no chance 
for amorphization to take place: the cost of breaking the 
strong ionic bonds is just too large. In fact, the system 
is found to undergo a transition to a complex — but or- 
dered — phase that maintains the chemical order of the 
system, i.e., that introduces no wrong bonds. Thus, it ap- 
pears that, indeed, strongly-ionic materials are not good 
candidates to pressure-induced amorphization: more "vi- 
olent" processes are required, and this suggests that the 
amorphous phase cannot exist in absence of wrong bonds. 

There is definite experimental evidence that MeV-ion- 
bombardfid c-InP contracts with respect to equilibrium 
material.Ej The density of a-InP, in fact, is slightly (a 
fraction of a %) larger than that of c-InP, which is a bit 
surprising in view of the fact that both a-Si and a-CiaAsj 
are less dense than their crystalline counterparts. Ld^Ell 
Our calculations are consistent with this observation in 
that the system is found to be, on average, overcoordi- 
nated, while the average bond length is larger than in the 
crystal. In contrast, using ab-initio MD and TB-MD, we 
have found a-GaAs — which is less dense than.|-c,-|G|aAs 
— to be undercoordinated, albeit only slightly.aEja 

Our paper is organized as follows. In Sec. |l|, in or- 
der to put our contribution in proper perspective, we 
present a summary of the information kno wn f rom both 
experiment and theory on a-InP. In Sec. Ill, we pro- 



vide details of our computational framework, including 
a description of the ground state (zinc-blende) and high- 
pressure (NaCl) phases. Discussion of our results for the 
model prepared by melt-and-quench is given in Sec. IV. 



There, we first present the results for the liquid phase, in 
particular static structure and diffusion. The structure 
and properties of the amorphous, low-temperature, phase 
is discussed next in terms of radial distribution functions, 
static structure factors, distribution of bond and dihe- 
dral angles, coordination numbers and bonding charac- 
teristics, vibrational properties, and density of electron 
states. In Sec. ^ finally, we present our findings on the 
possibility of amorphizing InP through the application of 
pressure. 



II. BACKGROUND 
A. Experiment 

To our knowledge, only very fey-X-ray or EXAFS,I 

and only one electron diffraction,ll3 experiments of a-InP 
(prepared by flash evaporation and deposition, and usu- 
ally non-stoichiometric) have been reported; the accuracy 
of these meaurements, as we discuss now, does not allow 
a precise determination of the local atomic order. The 
measured nearest-neighbour distances and coordination 
numbers for the various types of correlations are listed in 
Table ||. The error bar on the nearest-neighbour distances 
is reported to be ±0.02-0.05 A, while on the coordina- 
tion numbers, these are of the order of ±0.4-0.5, but this 
also depends on the model used to fit the EXAFS data, 
as can be seen in Table | for InaaPgv. 

The proportion of wrong bonds in the above measure- 
ments is reported to be anything between 10 and 40 % 
(after taking care of the off-stoichiometry of the sam- 
ples). The origin of the wrong bonds is not at all clear; 
part of the problem arises because of variations of com- 
position (which, p± this level, can be very significant): 
while Flank et alx3 believe that the system partly phase 
separates (i.p,-i clustering of the excess P takes place), 
Udron et alB2 indicate that the P is more-or-less uni- 
formly distributed in the samples. It is also suggested 
that wrong bonds are due to local composition fiu(jtii*u 
tions rather than the presence of topological defects ji3'E^ 
and in particular odd-membered rings. 

The total coordination number of each species is found 
in general to be quite close to 4, as can be seen in Ta- 
ble I (summing the partial coordinations), and rthis is 
also consistent with core-level-shift measurements.113 It is 
therefore tempting to conclude that odd-membered rings 
are absent in a-InP (though, of course, a perfect average 
total coordination of 4 does not preclude the existence 
of odd-membered rings) and, likewise, that wrong bonds 
are absent. However, using the same setjof EXAFS data 
but different fitting schemes. Flank et altB found the co- 
ordination of In to be either 4 or 4.8, as indicated in 
Table |. The error bar on these numbers is very signifi- 
cant and evidently prevents firm conclusions from being 
drawn. In fact, based on the "measured" coordination 
numbers. Flank et al. find a proportion of wrong-bonded 



2 



In atoms of almost 40%, but only 10% for P; taking the 
stoichiometry of the sample into account, they indicate 
that their fitting model must be incorrect. One must 
conclude, therefore, that the error bar on the experimen- 
tal values of the partial coordination numbers is so large 
that, for all practical purposes, they are at present un- 
known. 

Electron— I diffraction experiments have also been 
performedlZl on a-InP samples also prepared by flash 
evaporation and deposition. Though such quantities as 
partial coordination numbers arc not provided, a detailed 
analysis of the total pair correlation function of the mate- 
rial suggests that a-InP is "more disordered" than a-Ge 
and other III-V's: the first nearest-neighbour distance 
is larger that its crystalline counterpart, with a rather 
wide spread in the distances; the second peak is shifted 
towards smaller values, indicating that the average bond 
angle is smaller; finally, the third peak is more-or-less 
buried in the background, suggesting that order is to- 
tally lost beyond second nearest-neighbours. 

These electron diffraction measurements were inter- 
preted in terms of "standard" ball-and-stick, contiiipous- 
random-network models, namely those of PolkcS and 
Connell and Temkin.E3 The Polk model contains odd- 
membered rings, while the Connell- Temkin model does 
not. The measured structure of a-InP seems to be more 
adequately described by the unrelaxed Council- Temkin 
model, i.e., without odd-membered rings, in line with the 
above remarks. (This also agrees with a recent study of 
the structure of a-GaAs, as discussed in the next section). 
It is expected, however, that relaxation of the Connell- 
Temkin model would bring about odd-membered rings. 



B. Theory /Models 

To our knowledge, no structural model specific to a- 
InP has been ever been proposed. Only generic ball-and- 
stick models (Polk, Connell- Temkin) have been used to 
interpret strutural data; no computer model, based on 
any kind of potential, has been reported. 

Based on such a generic model, the density of electroD. 
states has been calculated by O'Reilly and Robertson.Efl 
It has been found that wrong bonds, presumably the 
most important type of defects in this material, lead to a 
significant number of states in the gap, and are therefore 
extremely important in determining the electronic prop- 
erties of the material. We will bring additional evidence 
for this in the present paper. 

It is appropriate to mention at this point that an opti- 
mized model for another III-V compound, a-GaAs, was 
very recently dpAcloped by one of the authors (LJL) and 
a collaborator .Eil LIsing the new "activation-relaxation 
technique" (ART)E3 for relaxing complex structures at 
K, a model was built that possesses almost perfect co- 
ordination and is essentially free of wrong bonds. In this 
study, it was demonstrated that a Connell- Temkin-like 



model, which contains no odd-membered rings, provides 
a better description of o-GaAs than a Polk-type model, 
which is more appropriate to elemental semiconductors. 
Thus, odd-membered rings must be present in a-Si and 
relatively rare in a-GaAs. How these conclusions apply 
to InP, however, is not clear: In addition to ionicity, InP 
differs from GaAs in that the atoms are significantly dif- 
ferent in size — Ga and As belong to the same row of 
the Periodic Table, while In and P are two rows apart. 
Thus, the competition between elastic deformation en- 
ergy and Coulomb repulsion will be rather different in 
the two materials. It is not possible, at present, to carry 
out ART simulations in order to address this issue since 
there exists no model potentials for InP. 



III. COMPUTATIONAL FRAMEWORK 

As mentioned earlier, our calculations were carried opt, 
using now standard first-principles molecular-dynamicsES 
plane- wave/pseudopotential methodology in the local- 
density approximation (LDA),lZI with the exchange-i 
correlation term expressed in the Ceperley- Alder form.ca 
The version of the code we use, however, is an imple- 
mentation of it optimized to run on a block of 32 nodes 
on a massively-parallel Cray T3D computer located at 
EPFL. As discussed below, this has allowed us to carry 
out extremely long runs in comparison to what would 
have been possible on a scalar machine. 

All calculations were performed on a constant-volume 
64-atom supercell for stoichiometric InP, i.e., 32 In and 
32 P atoms. The supercell volume was however changed 
"by hand" when appropriate (see below). The plane 
waves were cutoff in energy at 12 Ry, which proves to 
be essentially converged as far as structural properties 
are concerned, according to our tests (see also Rcf. |9|). 
The interaction between electrons and ion cores is de- 
scribed in terms of norm-conserving, fully-separable ah. 
initio pseudopotentials of the Kleinman-Bylander form.o 
Only the F point was used for integratine_the Brillouin 
zone. The program uses preconditioning so a rather 
large timestep of 10.0 a.u. (about 0.25 fs), with a cutoff 
"mass" of 3.0 a.u., could be used. The fictitious mass 
of the electrons was set to 300 a.u. A Nose thermostat, 
with a "mass" of 4.32 x 10^" a.u. was used to control the 
temperature; we have verified that the structural and 
dynamical properties of our systems are not significantly 
influenced by this choice. 

We have calculated the total energy as a function 
of lattice parameter using the above model for both 
the zinc-blende (ZB, F43m) and the sodium-chloride 
(NaCl, Fm3m) structures. The latter structure corre- 
sponds to thiLhigh-pressure phase of InP and other III-V 
compounds.OCJ The results are shown in Fig. |l|(a). The 
total-energy data were fitted to the "universal binding- 
energy function" (see, e.g., Ref. Eil) 
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where og is the fitted equilibrium lattice parameter and 
a and /3 are other fitting parameters. Fig. |^(b) shows the 
pressure (P = —dE/dV at K) for the two phases. 

From the data of Fig. |l|(a), we find the lattice param- 
eter to be 5.68 and 5.24 A for the ZB and NaCl struc- 
tures, respectively. For the ZB structure, the calculated 
lattice parameter is about 3% smaller than the experi- 
mental value (5.859 A). This discrepancy is largely due to 
our use of the LDA (which systematically underestimates 
lattice parameters), and to some extent also to limited 
Brillouin-zone sampling: the "trueH, LDA value, obtained 
by detailed integration, is 5.74 A,l3 2% smaller than ex- 
periment. For the NaCl structure, the computed value 
for the lattice constant using the same computational pa- 
rameters as above is about 9% smaller than experiment 
— 5.24 vs 5.71 A. Accurate bulk calculations on this 
phase show that the error originates in part (about 4%) 
from insufhcicnt Brillouin zone sampling (F-point only) 
and no Fermi-energy smearing scheme (the NaCl phase 
is found to be metallic at the theoretical equilibrium vol- 
ume). A further 4% of the error is recovered by using the 
non-lineat_jCore correction for the exchange-correlation 
potential,c3 leaving a residual error of about 1% due to 
the LDA and pseudopotential approximations. In the 
light of these results, we cannot expect our model to pro- 
vide an accurate description of this phase with the run- 
time calculation parameters reported above, which were 
required for the very long production simulations needed 
(see below). However, since we are primarily interested 
in the ZB phase, this will be of relatively little conse- 
quence, and we will still be able to draw qualitative con- 
clusions on the possibility of pressure-amorphizing InP 
(Sec. V). The energy difference between the two phases 
in our calculation is found to be about 0.11 eV/atom (in 
favor of ZB), compared to the fully-converged value of 
0.15 eV/atom and to about 0.38 eV/atom from experi- 
ment. While the error bar on the experimental value is 
not known, it is likely that part of this difference is due 
to the LDA approximation. 



IV. MELT-AND-QUENCH AMORPHIZATION 

A. Thermal cycle 

The thermal cycle used to prepare the amorphous sam- 
ple by melt-and-quench is summarized in Fig. |: Starting 
with a perfect crystal, the system was first equilibrated 
at room temperature (300 K), then heated up in steps 
until it melted, and finally cooled as slowly as possible 
into a glass. It should be stressed that the cooling rate 
used here, about 2 x 10^^ K/s, is probably the smallest 
ever achieved in a first-principles simulation of the liquid- 
glass transition: the total time covered is a formidable 90 



ps, compared to, typically, ~10 ps in corresponding sim- 
ulations of other materials. In spite of this, effects of 
the finite (and still large) cooling rate are expected to be 
present. 

For the lattice parameter of the crystalline phase, we 
used, at all temperatures, the value obtained above from 
the K global optimization; it should bejjiatpd that the 
thermal expansion of c-InP is very small,LjL3 and there- 
fore neglecting this effect is of little consequence. For 
the liquid, now, the density is larger than that of the 
crystal. This quantity is very difficult to calculate in 
the present simulation scheme, but is known (approxi- 
mately) from experiment. Thus we used, for the liquid, 
the density of the crystal scaled up by a factor equal ta 
the experimental ratio of liquid-to-crystalline densities,c£l 
namely ^^5. 1/4. 77=1. 069. The amorphous phase, finally, 
is known from experiment (on ion-implanted a-InP) to 
have a dcnsityj_almost exactly equal (to within 0.5%) to 
that of c-InP.o Again, here, this quantity is very dif- 
ficult to calculate in the absence of a constant-pressure 
option; thus, we simply assumed the amorphous-phase 
density to be the same as that of the crystal, an approx- 
imation which should be insignificant compared to other 
limitations of the study. 

Upon heating, the density of the system was changed 
from that of the crystal to that of the liquid at 1800 K, 
i.e., somewliai above the experimental melting tempera- 
ture of InPja Tra = 1335 ± 50 K. We found the system 
to remain crystalline at this temperature, i.e., to be in 
a super-heated state, a consequence of the finite (short) 
run time. It was then heated up to 2100 K, and found to 
melt, and then to 2400 and 3000 K, the highest temper- 
ature considered in this study. After coohng (in steps) 
to 2100 K, the density was changed back to that of the 
crystal, and the system "annealed" at 2400 K so as to 
remove the effects of the change in density. Quenching 
into the glass was then carried out in steps, proceeding 
more and more slowly into structural arrest (see Fig. |^.) 

The system was found to remain liquid (non-zero diffu- 
sion on the timescale of the simulations) at temperatures 
as low as 900 K, indicating a rather strong hysterisis of 
the melt-freeze cycle. While this is likely a manifestation 
of finite run times, it can also be attributed, in part, to 
our use of the LDA, which tends to underestimate the 
temperature of such transitions: For instance, in a frefe. 
energy calculation of the melting of Si, Sugino and CarE2l 
found a transition temperature somewhat below (300 K) 
that observed experimentally. It is however expected that 
finite-size effects on the transition temperature are rela- 
tively small. The liquid-glass transition can be seen very 
clearly in Fig. ^, which shows the total energy versus tem- 
perature upon going through the transition at constant 
density. 

The ground-state energy of the amorphous phase lies 
approximately 0.24 eV/atom above that of the crystal. 
This quantity (the latent heat of crystallization) has to 
our knowledge never been measured in InP; for Si, it 
varies between 0.14 and 0.20 eV/at, depending on the 
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state of relaxation of the material. Because of the pres- 
ence of wrong bonds, the heat of crystallization is ex- 
pected to be larger in III-V materials than in elemental 
semiconductors, consistent with our result. 



B. Liquid phase 

The liquid was studied in detail at four different tem- 
peratures: 3000, 2700, 2400, and 2100 K. The calculated 
diffusion constants are presented in Table ||. The er- 
ror bar on these numbers is estimated to be of the order 
of 10%, arising mainly from the limitations inherent to 
the method (size and time) . We find no significant differ- 
ences in the diffusional behaviour of the two components. 
From these data, we find an activation energy of about 
0.35 eV. To our knowledge, the diffusion constants are 
not known from experiment for InP; to giiie an experi- 
mental reference for a comparable system,c3 in the case 
of liquid GaAs, Z? = 1.6 x 10""* cmVs at 1550 K, i.e., a 
bit larger than the values we find here (extrapolating to 
lower temperatures). 

The structure of the liquid at the various temperatures 
considered was analyzed in terms of radial distribution 
functions, static structure factors, and coordination num- 
bers. The velocity auto-correlation functions and distri- 
bution of vibrational states were also calculated. 

The partial radial distribution functions (RDF's) 
gij{r) — pij{r)/ATrr'^CiPo (where Pij{r) is the correlation 
function for i-j pairs, is the relative concentration of 
species z, and po is the average number density) pro- 
vide detailed information about the short-range arrange- 
ments of atoms in the system; they are shown in Fig. 
^ for the lowest T examined at the liquid density, viz. 
2100 K. We find the hquid, independently of temper- 
ature, to have relatively little structure, essentially re- 
stricted to the first or perhaps second nearest-neighbour 
peak. Thus, there are essentially no correlations beyond 
a distance of about 3.5 A, and the "minimum after the 
first peak" is almost non-existent, except for P-P corre- 
lations, which seem to exhibit a well-defined minimum as 
well as a second-neighbour peak at this temperature. As 
discussed below, this absence of a marked structure will 
make it rather difficult to define coordination numbers. 

Likewise, we show in in Fig. ^ the partial, Sij{k), and 
total, S{k), static structure factors (SSF's) of the liq- 
uid at the lowest temperature. The SSF's are related to 
the RDF's by a Fourier transform and are in principle 
available directly from scattering experiments (neutrons, 
x-rays, etc). The SSF's were evaluated directly in re- 
ciprocal space in order to avoid the spurious oscillations 
that arise in the Fourier transform of a function that 
does not terminate smoothly (as is the case for finite- 
size models). Just like the radial distribution functions, 
the static structure factors show relatively little struc- 
ture. We know of no experimental Z-InP data to compare 
these results with. 



The results of Fig. ^ show the most strongly-marked 
correlation at short range to consist of In-P hetero- 
bonding. This is roughly twice as important as In- In and 
P-P bonding, which are nevertheless present in very sig- 
nificant number. Thus, "wrong bonds" are very present 
in this phase (and of course totally absent in the per- 
fect crystal), very likely a consequence of the metallic- 
bonding properties of the liquid, and in qualitative agree- 
ment with thfi first-principles calculations of Z-GaAs by 
Zhang et alM 

As noted above, defining coordination numbers in such 
a system is not simple. We plot in Fig. ^ the "running co- 
ordination numbers", i.e., integrated radial distribution 
fimctions, Zij(r) = pij(r)dr. If coordination numbers 
were well defined, one would see "plateaux" in these func- 
tions, corresponding to the successive neighbour shells, 
i.e., minima in the corresponding radial distribution func- 
tions. Clearly there are no such plateaux here. Never- 
theless, we list in Table HI the coordination numbers 



obtained by choosing some "reasonable" first-neighbour 
distances (as indicated in the Table). 

We find, despite the large error bars, the coordination 
numbers to decrease markedly with decreasing tempera- 
ture, i.e., the covalent character of the material is increas- 
ing upon approaching the transition temperature. This 
is true of all three types of partial correlations, and of 
course also of the average (total) coordination number. 
The latter, Z, decreases from 8.7 at 3000 K to about 7.0 
at 2100 K. We can extrapolate that, at the melting tem- 
perature of InP (1335 K), Z would be about 6.0, as is 
approximately found in Si just above melting. 

Fig. 0, finally, gives the density of vibrational states 
5(1^) for each atomic species, as well as overall. These 
were obtained by Fourier transforming the velocity auto- 
correlation functions. Though the density of states for 
In atoms show essentially no structure — it decreases 
rapidly with frequency — that for P atoms possesses a 
shoulder in the 20-40 meV range. This is likely related 
to the transverse and longitudinal optical phonon peaks 
in c-InP, respectively at 41 and 45 meV (in the present 
model; see below), and is manifest of the fast motion of 
the light P against the heavy In atoms. The frequency 
of this peak should therefore increase upon decreasing 
the temperature; indeed, this is is what we find upon 
examining g^iy) at various temperatures (not shown). 



C. Amorphous phase 

1. Radial distribution functions and static structure factors 

The partial RDF's for the fully-relaxed amorphous 
model at 300 K are presented in Fig. ||; also shown is the 
total (equi- weighted) RDF, g{r). We observe that the 
partial In-P RDF is quite similar (in shape) to the total 
RDF, reflecting the fact that, as expected, unlike-atom 
correlations largely dominate in the amorphous sample 
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at short distances. In the ideal ZB structure, of course, 
only hetero bonds are allowed and the first peak of the 
total RDF coincides with that of the In-P partial RDF. 
In the amorphous material, homo bonds are possible to 
some extent, even though hetero bonds prevail, as we 
discuss below. 

The presence of homo bonds is especially evident in 
the P-P partial correlation; they manifest themselves as 
a small peak in the RDF at a distance of 2.19 A, close 
to the P-P covalent bond distance (2.20 A, twice the co- 
valent radius). This distance is somewhat shorter than 
the In-P bond distance (2.51 A). In the case of In-In, we 
observe a shoulder, or prepeak, at a distance of 2.81 A, 
now larger than the In-P bond distance, but again close 
to the covalent bond distance (2.88 A). These effects can 
clearly be attributed to the size and ionicity differences 
between the two species (In is substantially larger than 
P.) In contrast, in a-GaAs, like-atom peaks are found at 
about the same distance as the unlike-atom peak. The 
nearest-neighbour distances we find agree quite closely 
with those from experiment reported in Table | (though 
at different chemical compositions). For P-P, we find 2.19 
A, vs 2.20-2.24 from experiment; for In-P, we obtain 2.51 
A, compared to 2.57-2.59 experimentally; and for In-In, 
which is most difficult to define, as is also the case experi- 
mentally, we have 2.81 A vs 2.76-2.98. We note that part 
of the observed difference arises from our model underes- 
timating (by about 3%) the lattice parameter of the real 
material as discussed earlier; in view of this, we conclude 
that our model is in close agreement with experiment 
as far as nearest-neighbour distances are concerned and 
modulo the error bars inherent to both methods. 

In the crystal, the equilibrium LDA In-P bond dis- 
tance is 2.46 A, while second nearest-neighbours lie at 
4.02 A. In our amorphous sample, we find, from the to- 
tal RDF (Fig. ^), the nearest-neighbour peak at 2.51 A, 
a bit larger than the corresponding value in the crystal. 
In contrast, the second-neighbour peak is at about 3.9 
A, thus shifted towards smaller values compared to the 
crystal, and is much broader. In fact, it is clear from 
Fig. ^ that the second peak is made up of at least two 
subpeaks, with a shoulder at about 4.4 A arising from 
In-P correlations. In any case, these observations agree 
with the electron diffraction data of Ref. |l^, discussed in 



Sec. [I A 



It is clear from Fig. g that the concept of nearest- 
neighbour distances in the amorphous phase is somewhat 
ill-defined, especially in the case of In-In correlations, 
where the first peak is almost merged into the second 
peak. In fact, the "second- neighbour" peak, for all cor- 
relations, is rather wide, consisting of several sub-peaks, 
refiecting the large spectrum of possible configurations in 
the disordered phase. All correlations seem to differ little 
from unity beyond the second peak, indicating that or- 
der, in the amorphous phase, is indeed very short-range, 
restricted to the first- and second-, perhaps third-, neigh- 
bour shells. We note, also, that second-neighbour peaks 
differ very significantly in shape from the corresponding 



peaks in crystalline material. 

For completeness, we present in Fig. ^ the partial 
and total SSF's for our model sample. The total SSF 
was obtained by combining the partial Sij{k) with equal 
weights. [In principle, >S'(fc) is a weighted sum of the 
partials, where the weights are related to the scatter- 
ing lengths of the atoms for the probe used.] The total 
interference function (essentially the SSF) of a-InP at al- 
most stoichiometric concentration, measured by electron 
diffraction, has been reported by Gheorghiu et al.tll They 
observe a small peak at 2.1 A~^, and three large peaks 
at 3.5, 5.7, and 8.0 A~^, respectively. This correlates 
extremely well with the total SSF displayed in Fig. ^. 



2. Bond and dihedral angles 

We give in Fig. ^(a) the distribution of bond angles 
in the amorphous structure, all combinations taken into 
account. As can be inferred from the above discussion, 
the definition of "bond" is somewhat arbitrary. The cut- 
off distances we used, extracted from the corresponding 
RDF's (Fig. |) are 2.91, 3.13, and 2.55 A for In-In, In-P 
and P-P, respectively; the value for In-In, which hardly 
exhibits a nearest-neighbour peak, is subject to a signif- 
icant error. These cutoff distances will also be used for 
determining the coordination numbers, below. 

The bond-angle distribution is wide but exhibits a 
strong peak at about 107°, a little bit smaller than the 
tetrahedral angle (109.5°). A shift of the bond-angle 
peak to smaller values has also been observed by electron 
diffractionO A similar shift has been obtained theoreti- 
cally by Stich et al. for a-Si.l3 The bond-angle distribu- 
tion here differs from the case of a-Si in that it shows 
a rather marked shoulder at about 90° — likely arising 
from four-membered rings and from those atoms that are 
five- or six-fold coordinated — as well as a weak shoulder 
at about 150°, which perhaps originates from three- fold 
coordinated atoms. It is quite remarkable that there ex- 
ists almost no correlations with an angle of 60°. This 
is in sharp contrast with other tetrahedral semiconduc- 
tors (elemental or compound) , modeled either empirically 
or using TB or first-principles MD, where a significant 
peak or shoulder is observed at such small angles, aris- 
ing from small, e.g., three-membered, rings. This indi- 
cates that the chemistry of this system is robust enough 
that such defects are rare (three-membered rings are ex- 
tremely costly in both elastic-deformation and electronic- 
repulsion energies while four-membered rings cost only 
elastic energy), as can indeed be verified in Table IV, 
and/or the relaxation of the present model has been par- 
ticularly effective. (Ring statistics are extremely sensi- 
tive to details of the local structure, and in particular co- 
ordination; this explains the sizable differences between 
a-InP and a-GaAs in Table |rv|.) 

In Fig. |l^(b), we give the distribution of dihedral an- 
gles (angles between second- neighbour bonds). In the ZB 



6 



structure at low temperatures, the corresponding distri- 
bution has two sharp peaks, at 60 and 180°. In the case 
of a-InP, we observe a rather flat distribution, except for 
two small dips at and 120° (which are equivalent, on av- 
erage, for tetrahedral systems), perhaps a memory of the 
crystalline phase, but in any case much less pronounced 
than the corresponding ones in a-GaAs,t3 which chemi- 
cally orders a bit more strongly than a-InP (see below). 



part, for the observation of a lower density in a-GaAs 
than in c-GaAs.Ej (Disorder itself is expected to cause 
a decrease in density). In contrast, the predominance of 
overcoordinated defects in a-InP is likely responsible for 
its larger density compared to c-InP,E3 given, as we have 
seen above, that the average bond length in the amor- 
phous phase is larger than in the crystal. 



3. Coordination numbers 

The average coordination numbers can be obtained by 
integrating the appropriate RDF's up to the nearest- 
neighbour distances defined above; this and other rele- 
vant numbers are listed in Table while the "running" 
coordination numbers — now exhibiting plateaux — are 
presented in Fig. |Tl|. We obtain in this way a total co- 
ordination number of Z = 4.27, in reasonable agreement 
with the experimental value mentioned above, Z « 4, 
i.e., within the uncertainties inherent to both methods. 

A detailed picture of the short-range structure is pro- 
vided by the partial coordination numbers Zi and Zij, 
I, J = In or P, also hsted in Table 0. We find the partial 
coordination numbers of In and P to be almost identi- 
cal — 4.25 and 4.29, respectively. Modulo the limita- 
tions mentioned above, this is again in agreement with 
the available experimental values (4.0 < Zi^-in < 4.8 
and Zp.p w 4.0). Thus, despite the large difference in 
size, and because of the strongly ionic character of the 
material, each atom is surrounded by the same average 
number of atoms. In this sense, it can be said that all 
atoms occupy the same volume. 

If we detail further the average coordination numbers, 
we find, from Table ^ that coordination essentially con- 
sists of hetero bonding, i.e., the system is chemically or- 
dered. Thus, in the case of In, out of the 4.25 neighbours, 
3.91 are P and only 0.34 are In. Likewise, for P, which 
has 4.29 neighbours, we have 3.91 In and 0.38 P. We also 
see, upon comparing with the TB-MD results for a-GaAs, 
that the chemical short-range order appears to be a bit 
stronger in the latter. However, it must be said that the 
data reported in Table ^ (and following) were obtained 
using the ART procedure, which allows more extensive 
relaxation of the network than is possible with MD. 

Modulo the error bars, the average coordination num- 
ber is larger in the amorphous phase than in the ZB crys- 
tal at equal density, i.e., there are a number of overco- 
ordinated atoms. This can be seen in Table VI, where 



we present the distributions of coordination numbers in 
our amorphous sample. Even though the distribution is 
rather sharply peaked, there are nevertheless a signifi- 
cant number of coordination defects. In fact, we find, 
overall, very few (1.9%) atoms that are undercoordinated 
{Z < A), while quite many (26.8%) are overcoordinated. 
This, again, contrasts quite sharply with o-GaAs, which 
is slightly undercoordinated; this might be the cause, in 



4-. Chemical disorder and wrong bonds 

A quantitative measure of chemical correlations in the 
binary compound AB is provided by the "concentration- 
concentration" coordination number. Zee — cb{Zaa ~ 
Zba) + ca{Zbb - Zab) (see, e.g., Ref. ||) where Ci is the 
concentration of i-type atoms in the system. Zee = —4 
exactly in c-InP; for our amorphous sample, we find 
Zee = —3.55 (Table indicating, as was already ev- 
ident from the above discussion, a rather strong chem- 
ical order. Chemical order can also be quantified in 
terms of the generalized Warren chemical short-range 
order parameterJ3 aw = Zeej {cbZ a + caZb^, where 
Zi = 'Ylij . aw = indicates complete randomness 
whereas positive and negative values indicate preference 
for homo and hetero nearest-neighbour coordination re- 
spectively. Evidently, in c-InP, aw = —1; for a-InP, 
we obtain aw = —0.84 (cf. Table 0), revealing, again, 
a strong preference for chemical ordering, a bit weaker, 
perhaps, than in a-GaAs, for which the TB-MD model 
gives Zee = —3.54 and aw = —0.88. This, again, reveals 
the importance of Coulombic ordering in InP and GaAs. 

The overall similarity between the RDF's of group- IV 
materials and the III-V semiconductors suggests jfeat the 
materials have comparable short-range structured How- 
ever, as discussed above, there exists a significant number 
of coordination defects, such that the overall coordina- 
tion exceeds, in the present case, the canonical value of 
4. Likewise, the structure exhibits a significant number 
of "anomalous" rings — as can be seen from Table IV 
— and in particular odd-membered, just as they can be 
found in a-Si or a-Ge. An immediate consequence of this 
is that there must exist "wrong" bonds in the structure. 
We find in our model that 8.4% of the bonds arc wrong 
(cf. Table Such a proportion of wrong bonds is re- 
markably small in view of the fact that the system is 
slightly overcoordinated and thus is manifest of the ex- 
cellent quality of the model. 

Experimentally, the proportion of wrong-, bonds has 
been reported to lie in the range 10-40%.t3 The large 
spread in the values reported is explained by the fact that 
some samples are believed to phase separate. It has been 
conjectured, also, that the wrong bonds in a-InP might 
originate from local composition fluctuations rather than 
coordination defects. Our calculations indicate that co- 
ordination defects are responsible for the wrong bonds. It 
is perhaps appropriate to remark that it is quite difficult 
to imagine an amorphous network without coordination 
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defects and/or odd-membered rings, but of course the 
density of such defects is not known precisely and prob- 
ably depends quite strongly on the "method of prepa- 
ration", be it experimental or computational. In fact, 
the "effort" required to reduce the proportion of wrong 
bonds to a value smaller than the present 8.4% would 
appear to be formidable if it is a consequence of model 
limitations, and in particular the quench rate used in the 
MD simulations. In any event, the proportion of wrong 
bonds we obtain here must be taken as an upper limit to 
the actual value: it is certainly the case that the num- 
ber would decrease if corresponding simulations were car- 
ried out on a larger system, so as to minimize the elastic 
constraints, on longer timescales, in order to allow more 
complete Mlaxation. An ART optimization could resolve 
the issueJlil'LS however, this is presently not feasible ab 
initio^ or otherwise, since there exists no model potentials 
for InP. 



one, In. (Cf. also the discussion concerning Fig. ^ in Sec. 
IVB .) In view of this, and of the fact that wrong bonds 
do not exist in elemental semiconductors (while other de- 
fects, e.g., coordination, do), we conjecture that the band 
at 55 meV arises from the motion of phosphorus atoms 
against one another, i.e., P-P wrong bonds. (Because of 
the heavier mass of indium, In-In wrong bonds will show 
up at much smaller energies, and therefore be buried in 
the continuum of states. For similar reasons, a high- 
energy wrong-bond peak has not been observed in the 
DOS of a-GaAs.)liJ It would be of utmost interest that 
experimental confirmation of this point be carried out, 
since this would give a direct indication of the presence 
of wrong bonds and a measure of their relative impor- 
tance. 



6. Electronic properties 



5. Vibrational properties 

The partial and total densities of vibrational states 
(DOS) as deduced from our model are presented in Fig. 
0(a); for reference, we give, in Fig. |l|(b), the corre- 
sponding DOS for c-InP calculated within the same com- 
putational framework. To our knowledge, there exists no 
experimental measurements of this quantity for a-InP, 
while the vibrational sppc.t.rum of the crystalline mate- 
rial is well characterized. OtZI Thus, in c-InP, rather wide 
TA and LA bands are found in the range 6-9 and 20-23 
meV, respectively, while more sharply defined, TO and 
LO peaks are seen at about 37 and 42 meV, respectively. 
In the present calculation, we find a large TA peak at 
about 7-8 meV, and a fairly broad LA band in the range 
16-22 meV. The sharp optic peaks are found at 41 and 45 
meV, respectively. Thus it appears that our model over- 
estimates a little bit the energy of the optic peaks, while 
it underestimates a little the energy of the LA band. It 
must be said, however, that the low-frequency acoustic 
modes are the most difficult to probe with molecular dy- 
namics (explaining, in part, the oscillatory structure at 
low energies). 

The density of states of our model a-InP agrees, 
"broadly" speaking, with that of c-InP, except for a sig- 
nificant softening of the higher-energy peaks. The total 
DOS exhibits a well-defined peak at about 7-10 meV, 
corresponding to the crystal's TA peak, a broad band 
centered at about 18 meV, close to the crystal's LA peak, 
and two well-defined peaks at 32 and 38 meV, corre- 
sponding to the crystal's TO and LO bands. 

Fig. |l2| reveals yet another feature in the DOS which 
is absent in the crystal, as well as in elemental semicon- 
ductors, namely a (rather broad) peak at high frequency 
— about 55 meV. It is clear from Fig. |l^ that the optic 
peaks are primarily associated with the fast and ener- 
getic motion of the lightest atom, P, against the heavier 



We have calculated the density of electron states g{E) 
for our model a-InP at K; the results are shown in Fig. 
p^ . (In order to improve presentation — in view of the 
limited statistics of the model — the density of states has 
been smoothed lightly using a Gaussian filter of width 
0.15 eV.) The forbidden energy gap, about 1.08 eV, is 
clearly visible about the Fermi energy. For crystalline 
InP in the ZB phase, our computed value for the LDA 
direct gap at the T point is 1.50 eV (at a— 5.68 A) and 
compares well with the value 1.50 eV of Ref. ^ and with 
experiment, 1.42 eV. The gap of a-InP, therefore, is a bit 
smaller than that of the crystalline material. 

One important difference, however, is that there are 
defect states in the gap of a-InP not present in (ideal) 
c-InP. Thus, we have identified one particular electron 
level giving rise to contributions near mid-gap in g{E), 
clearly visible in Fig. |l3|. We have examined the local 
density of states for this particular level, and found that it 
corresponds to an empty, distorted, octahedral, mostly- 
indium (five out of six corners), "cage", i.e., basically, 
a cluster of wrong In-In bonds. That wrong bonds give 
rise to states in the gap has also been inferred from a 
comparison of a Polk-type (with wrong bonds) with a 
Connell-Temkin-type (without wrong bonds) model for 
o-GaAs, as discussed above (Sec. |II B| ). 

V. CAN AMORPHIZATION BE 
PRESSURE-INDUCED? 

In an attempt to verify the possibility that InP could 
amorphize under compression, we subjected the equilib- 
rium ZB structure to pressure by increasing "slowly" , 
in steps, the density. Referring to Fig. ||(b), we find the 
correspondence between pressure and density (in fact the 
lattice parameter). Starting with InP in its perfect ZB 
arrangement, properly equilibrated at 300 K, the lattice 
parameter was thus decreased (always at 300 K) from 
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5.68 A (equilibrium) to 5.23 A, i.e., down to a value 
smaller than (of. Sec. Ill) the equilibrium lattice constant 
for the NaCl phase. From Fig.|](a), we would expect the 
ZB crystal to undergo a transition to the high-pressure 
NaCl phase in the MD run when a is set to values be- 
low about 5.4 A. (This value corresponds to a pressure of 
about 140 kbar; at the highest density investigated here, 
for a — 5.23 A, the pressure in the ZB structure is 270 
kbars [cf. Fig. 0(b)].) 

Fig. ^ shows the evolution of the total energy of the 
system at 300 K relative to the ZB crystal at K, as a 
function of time-density. For reference, we also show on 
this plot the total energy of the equilibrium ZB phase at 
300 K, the total energy of the amorphous phase (also at 
300 K) at the equilibrium density, and the total energy 
of the NaCl phase at K, obtained as discussed earlier. 
We find that the energy increases rather smoothly with 
density. At a value of a = 5.33 A, visual inspection of the 
system indicates that it undergoes some sort of distortion 
into a state which is definitely not ZB, but which bears 
strong resemblence to it. This distortion is also visible 
in Fig. |lj as a slight decrease of the total energy as a 
function of time. At density values corresponding to lat- 
tice parameters smaller than about 5.28 A, the energy of 
the compressed ZB crystal exceeds that of the amorphous 
phase. Yet, no transition to an amorphous phase takes 
place. Upon increasing the density further, we observe 
another transformation for a = 5.23 A, clearly visible 
m Fig. |l| — the energy drops significantly, to a state 
which is evidently distinct from the NaCl structure (its 
computed energy is much higher). We have not analyzed 
this phase in detail but it is evidently ordered, as can be 
seen in Fig. and might possibly be an intermediate 
state on the way to the NaCl phase. One thing is clear, 
however: the new phase maintains the chemical order of 
the system, i.e., introduces no wrong bonds. It would 
seem, therefore, that pressure is not a proper route for 
amorphization; rather, the system prefers to reorganize 
into a new crystalline form, which is more favorable in 
view of the high cost in energy of wrong bonds. 

The argument might be presented in another way: In 
Fig. |l^, we see that the energy of the new phase, 0.33 
eV/at, lies only slightly above that of the amorphous 
phase obtained by the melt-and-quench cycle, 0.29 eV/at. 
The latter value, however, is at the equilibrium density. 
Under compression, the energy of the amorphous phase 
would also go up, presumably by an energy smaller than 
but comparable to 0.33 eV/at, and thus would largely 
exceed that of the compressed crystal. In view of this, 
it is very difficult to imagine that amorphization could 
take place under compression. Rather, a transition to 
the NaCl phase would take place. 

The emerging picture suggested by our simulations is 
that InP does not amorphize under pressure, even at val- 
ues largely exceeding those required for the system to 
transform into the high-pressure NaCl phase. The en- 
ergy of the amorphous phase lies well above that of the 
compressed crystal and the cost of breaking the strong 



ionic bonds is just too large. The system, rather, finds 
its way into a complex, but ordered, phase that main- 
tains the chemical order of the system, i.e., that does 
not introduce wrong bonds. It appears, therefore, that, 
indeed, strongly-ionic materials ar&Jiot good candidates 
to pressure-induced amorphizationtj more "violent" pro- 
cesses (such as implantation) are required, and this sug- 
gests that the amorphous phase cannot exist in absence 
of wrong bonds. Though we have not examined this, 
it is not impossible that InP would amorphize from the 
high- pressure phase upon the release of pressure and/or 
through proper heat/pressure treatment. In view of the 
above energy-wise arguments, however, this seems to be 
very unlikely. 

VI. CONCLUDING REMARKS 

We have carried a detailed and extensive first- 
principles molecular-dynamics study of the structure and 
electronic properties of amorphous InP produced by 
rapid quenching from the liquid. The structure of the 
material is found to be strongly ordered chemically, about 
the same, for instance, as in a-CaAs, even though there 
are a significant number of coordination defects (anti- 
sites) and despite the presence of odd-membered rings. 
We find, as a consequence, that there exists "wrong 
bonds" in the system, in an amount of about 8%; these 
are a consequence of the presence of defects in the system, 
not of composition fluctuations, as has been conjectured. 
The system, in fact, is found to be over-coordinated, 
which might be the reason for the observed higher density 
of a-InP compared to c-InP. We have also investigated 
the possibility of pressure-amorphizing InP. Our calcu- 
lations indicate that the cost of a transformation of the 
compressed zinc-blende crystal into an amorphous phase 
is so large that it is very unlikely that it would take place. 
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TABLE I. Experimental values of the nearest-neighbour 
distances and partial coordination numbers. For InsaPg?, the 
results of two different fits to the same EXAFS data are indi- 
cated. 



Sample 




rjvjv 


(A) 






Z 






Ref. 




In-In 


In-P 


P-In 


P-P 


In-In 


In-P 


P-In 


P-P 






a-InssPer 


2.98 


2.59 


2.58 


2.20 


1.5 


2.5 


1.8 


2.2 


















1.8 


3.0 








Li 




a-InssPes 


2.80 


2.58 


2.58 


2.24 


1.2 


2.8 


1.6 


2.5 








a-In4oP60 


2.76 


2.57 


2.58 


2.24 


0.9 


3.1 


2.2 


1.8 

































TABLE II. Diffusion constants in the liquid at various 
temperatures, in units of 10"* cm^/s. 



r(K) 


D 


3000 


2.13 


2700 


2.32 


2400 


1.96 


2100 


1.43 



TABLE III. Partial and total coordination numbers in the 
liquid; the cutofi' distances rz are also given. The total Z is 
obtained from the partials e& Z = . CiZij. 



T{K) 


In-In 


In-P 




P-P 




Total 




rz 


Z 


rz 


Z 


rz 


Z 


Z 


3000 


4.0 


5.6 


3.7 


5.1 


2.8 


1.6 


8.7 


2700 


4.0 


5.6 


3.6 


4.8 


2.8 


1.8 


8.5 


2400 


3.7 


4.2 


3.6 


4.8 


2.7 


1.5 


7.7 


2100 


3.5 


3.4 


3.5 


4.6 


2.7 


1.4 


7.0 



TABLE V. Structural properties of a-InP at 300 K: 
coordination numbers Z (partial, species, total, concentra- 
tion-concentration), Warren chemical short-range order pa- 
rameter OLw , and proportion of wrong bonds (WB). Also 
shown, for comparison, are the results for a-GaAs obtained 
from a fully-relaxed, ART-optimized, TB-MD modeLN Here 
A represents either In or Ga and B represents either P or As. 





Zaa 


Zab 


Zbb 


Za 


Zb 


Z 


"Zee 


aw 


WB 


c-InP 





4 





4 


4 


4 


-4.00 


-1.0 





a-InP 


0.34 


3.91 


0.38 


4.25 


4.29 


4.27 


-3.55 


-0.84 


8.4 


a-GaAs 


0.22 


3.75 


0.21 


3.97 


3.96 


3.96 


-3.54 


-0.88 


5.2 



TABLE VI. Distribution (in %) of total coordination num- 
bers for the a-InP sample at 300 K, as well as for the ideal 
c-InP structure. Also shown, for comparison, are the results 
for a-GaAs obtewned from a fully-relaxed, ART-optimized, 
TB-MD model.y 



z 





1 


2 


3 


4 


5 


6 


7 


c-InP 














100 











a-InP 











1.9 


71.3 


24.7 


2.1 





a-GaAs 











11.1 


82.8 


5.2 


0.6 


0.2 



TABLE IV. Number per atom of n-membered rings for the 
a-InP sample at 300 K, as well as for the ideal c-InP structure. 
Also shown, for comparison, are the results for a-GaAs 0|b, 
tained from a fully-relaxed, ART-optimized, TB-MD modelo 



n 


3 


4 


5 


6 


7 


c-InP 











4 





a-InP 


0.02 


0.44 


0.37 


2.35 


4.37 


a-GaAs 


0.05 


0.10 


0.21 


1.37 


0.76 
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FIG. 3. Total energy of the system through the liq- 
uid-glass transition; here, the density is that of the crystal. 
Also shown is the energy of the ZB crystal. The lines are for 
guiding the eye. 



FIG. 1. (a) Total energy and (b) pressure vs lattice par- 

ramctcr for the two structures considered: ZB and NaCl. The 
lines are obtained by fitting to Eq. (1) 
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FIG. 2. Thermal cycle of the melt-and-quench process 
used to prepare the amorphous phase, as discussed in the 
text. 
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FIG. 4. Partial and total radial distribution functions of 
liquid InP at 2100 K. For clarity, in this and similar figures, 
the zeroes are displaced along the y axis. 
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FIG. 5. Partial and total static structure factors of liquid FIG. 7. Densities of vibrational states of liquid InP at 

InP at 2100 K. 2100 K. 
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FIG. 9. Partial and total structure factors of a-InP at 300 



K. 



FIG. 11. Running coordination numbers corresponding to 
the radial distribution functions of Fig. H. 
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FIG. 10. Distribution of (a) bond and (b) dihedral angles 
in a-InP at 300 K. 
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FIG. 12. Densities of vibrational states of (a) a-InP and 
(b) c-InP at 300 K. 
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FIG. 13. Density of electron states for the amorphous 
sample at K. 
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FIG. 15. Top; Ball-and-stick representation of the final, 
compressed (a = 5.23 A), ZB crystal; bottom: same, with 
bonds removed in order to show better the underlying struc- 
ture. 
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FIG. 14. Energy versus time-density for the compressed 
ZB phase at 300 K relative to the ZB crystal at K. Also 
indicated are the energies of the NaCl phase at K, of the 
equilibrium ZB phase at 300 K, and of the amorphous phase 
at the equilibrium density. 
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